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Abstract: On contrary to the customary thought, the well-known "lemma" 
that the distribution function of a collisionless Boltzmann gas keeps invari- 
ant along a molecule's path represents not the strength but the weakness of 
the standard theory. One of its consequences states that the velocity dis- 
tribution at any point is a condensed "image" of all, complex and even dis- 
continuous, structures of the entire spatial space. Admitting the inability to 
describe the entire space with a microscopic quantity, this paper introduces 
a new type of distribution function, called the solid-angle-average distribu- 
tion function. With help of the new distribution function, the dynamical 
behavior of collisionless Boltzmann gas is formulated in terms of a set of 
integrals defined by molecular paths. In the new formalism, not only that 
the difficulties associated with the standard theory are surmounted but also 
that some of practical gases become calculable in terms of today's computer. 

PACS number: 51.10.+y. 
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1 Introduction 



The Boltzmann gas, by which we mean a gas in that molecules are subject 
to binary short-range interparticle forces, is usually studied by means of 
the Boltzmann equation Q [^. The conventional wisdom assumes that if a 
sufficiently large supercomputer were available, a finite-difference scheme 
could be employed to construct solutions of the Boltzmann equation and 
many realistic gases would become completely analyzable. 

Some of difficulties related to the solution-construction scheme of the 
Boltzmann equation are relatively well-known Q. They include the follow- 
ing, (i) Seven variables have to be dealt with (time, spatial coordinates and 
velocity components), (ii) The collisional operator in the Boltzmann equa- 
tion is complex in nature. Due to these difficulties, only low-dimensional 
problems, for instance those having one dimension of coordinates and one 
dimension of velocity components, were tried in practical calculations. 

In some of our works J§]|5|], it was argued that more fundamental diffi- 
culties existed with the standard kinetic theory. One of them was that dis- 
tribution functions of Boltzmann gases involved in discontinuity and quasi- 
discontinuity, thus the effectiveness of using a finite-difference scheme to 
solve the Boltzmann equation, as well as the validity of the Boltzmann 
equation itself, became questionable. 

Putting aside whether the arguments presented by our previous works 
are sufficiently convincing or not, we will, in this paper, engage our primary 
attention with the following subject. Is there any alternative approach with 
the aid of that one can use today's computational means to calculate at least 
some of full-dimensional practical gases? A positive answer to this question 
will apparently arouse widespread interest. 

It turns out that if an average-type distribution function, called the solid- 
angle-average distribution function, is introduced and the path-information 
of molecules is adequately taken into account, then the objective mentioned 
above can be accomplished. 

The structure of this paper is the following. In Sec. 2, the standard 
approach to the Boltzmann gas is employed to "derive" a lemma that the 
distribution function of a collisionless Boltzmann gas, if continuous at the 
beginning, will keep its value along a molecule's path. In applying the 
lemma to various situations, we are led to an unexpected conclusion that 
strictly continuous and well-behaved distribution functions exist only on a 
theoretical assumption. In Sec 3, we examine two kinds of discontinuous dis- 
tribution functions. The first kind is related to point-like molecular sources 
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and the second kind to surface-like molecular sources. It is shown that the 
solid-angle-average distribution function needs to be introduced in order to 
treat discontinuous distribution functions related to the molecular sources. 
Sec. 4 offers a comprehensive formulation in which both continuous and dis- 
continuous distribution functions are formulated in terms of path-integrals 
defined by initial-state and boundary conditions. The final section, Sec. 5, 
gives a brief survey on essential features of the new formalism. 



2 Evolution of continuous distribution function 

We suppose in this section that the distribution function of interest is per- 
fectly continuous, differentiable in a more exact language, in the position- 
velocity space (//-space). 

Under this supposition, according to the standard kinetic theory, the 
collisionless Boltzmann equation holds and it reads 

df df F df 
dt or m ov 

in which the collisional effects between molecules have been disregarded (as 
the title of this paper has suggested) . If the initial state 

/(t + 0,r,v) =/(t ,r,v) (2) 

and the boundary condition 

fit, r, v) = J K(y, vi)/(t, r, Vl )dvi, (3) 

where the shape of K represents the collisional effects between molecules 
and boundaries ||, are regarded as known ones, the solution of Eq. (||) can 
be constructed by means of a finite-difference schemeQ. In such a scheme, 
notably, the governing equation (|I|) and the initial and boundary conditions 
constitute a complete equation set, and any additional information, 
such as the path-information of individual molecules, becomes theoretically 
dispensable. (In a sense, the whole idea behind introducing the conservation- 
type partial differential equations into kinetic theory is just to avoid knowing 
paths of individual molecules.) 

Instead of making a detailed analysis of the standard formalism out- 
lined above, which has been done elsewhere we now use the formalism to 
"derive" a relatively well-known lemma in statistical mechanics. The point 
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here is that the lemma will automatically lead us to the very essence of gas 
dynamics. 

By relating the vector variables r and v in (jl]) to a molecule's path, the 
partial differential equation can be rewritten in the form of a simple ordinary 
differential equation 

I = »■ (4) 

at r(t),v(t) 

Of course, r(t) and v(i) are associated with the equations of motion for a 
single molecule 

dr(t) = dv(t) = F 
dt ' dt m K ' 

where F stands for all forces acting upon the molecule. Equation (Q) simply 
means 

/(*,r,v)| r ( t ) jV ( t ) = Constant. (6) 

That is to say, in the absence of collisions a continuous distribution function 
keeps invariant along a molecular path in //-space. [Since formula (|6|) is 
the integral of (||) along a path, we may regard the formula as a primitive 
prototype of path-integral approach.] If the time development operator T is 
employed to represent the equations of motion (||), namely if we have 

r(to,t)[r(«d),v(*o)] = [r(t),v(t)], (7) 

then we can equivalently express (||) as 

f(t, r, v) = f(t , tq, v ) = f(t ,T-\ T-V), (8) 

where ro = r(io), vo = v(io) and T _1 = T(t,to)- In view of (§) and (||), we 
will refer to /(to) as the "source" function and refer to f(t) as the "image" 
function. In a similar spirit, it will be said that ro is a source point and r 
an image point. 

Now, we use the path-invariance expressed by ([]) to derive some dynam- 
ical features of the Boltzmann gas. These features, though quite interesting 
and very essential, received no enough attention before. For simplicity only 
situations in that no external force exists will be considered, though equa- 
tions (|])-(|8]) are valid more generally. 

Suppose, to begin with, that the initial distribution function f(to,r,v) 
is known (which is continuous) and we are concerned with the velocity dis- 
tribution at a specific point r, denoted by / r (t,v) hereafter. Equation (||) 
shows that 

/ r (t,v) = /(t ,T- 1 r,v). (9) 
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For one fixed value of Vi = |vj|, the source points rj = T _1 r simply form 
a spherical surface, labeled as Si = S(vi) in Fig. 1. The radius of such a 
surface can be expressed by 

|r< — r| = Vi -(t-t ). (10) 

To determine the velocity distribution, we must let Vi take values from very 
small to very large, and thus the radius |r — r$| varies from very small to very 
large. This implies that the velocity distribution f r (t,v) is a condensed im- 
age of the entire spatial space. If f(to, r, v) is not spatially uniform, f r (t, v) 
will possess a complex structure, infinitely complex in general situations 
(which reminds us of the famous Cantor set and many other interesting 
structures in the fractal studies). To make the situation even more wor- 
risome, the inference presented above also suggests that any discontinuous 
things, if exist in the spatial space, will constantly create discontinuous im- 
ages on velocity distributions in the nearby and distant regions (virtually 
everywhere) . 

To see the point in another perspective, let's assume that a numerical 
work is being carried out and the velocity space is being divided into many 
small but finite cells (for time being pay less attention to the spatial space), 
of which one takes the form Av = v 2 Av AQ. It is then up to the investigator 
to determine the distribution function for each of the cells, denoted by 

f T (t,Av,AQ). (11) 

In other words, each of f r (t, Av,AQ) should be endowed with one definite 
value. The task, though seems quite straightforward, cannot be done con- 
ventionally. Look at Fig. 2, in which the spatial cone — AQ corresponding 
to the velocity cone AQ, called the effective cone hereafter, and the two 
spherical surfaces Si = S(vi) and SVk = S(vi + i), where vi + \ — v j = Av, 
are plotted. It is easy to find that when v varies within the range v 2 Av AO,, 
the source point T _1 r runs over the shaded volume in the figure. Since the 
thickness of the shaded volume is proportional to Av solely (for the fixed 
t — to), if we let Av be smaller and smaller the value of 

lim f r (t,Av,Q) (12) 
Av— >0 

tends to a limit in the usual sense. Whereas, the size of ASi, the area 
element on S(vi) enclosed by the effective cone, is proportional to vf, and if 
we let AO, be smaller and smaller the value of 

lim o f r (t,v, An) (13) 
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tends to a limit at a rate that strongly depends on the magnitude of v. 
That is to say, the velocity distribution at r is not uniformly continuous. 
According to mathematics textbooks Q [Q], a function that is not uniformly 
continuous may involve some irregularities and has to be treated carefully. 
As far as our situation is concerned, there indeed exist difficult things. As 
one thing, when V{ is rather large, the area of ASi must be rather large. 
If the source function /(to) ° n this large area element varies significantly 
(quite possible in practical situations), giving a unique value to the image 
function f r (t, Av, Af2) becomes a tricky business. This type of difficulty has 
been named as quasi- discontinuity in our previous works 

We now study boundary effects. At first suppose, if the path-invariance 
expressed by (||) is combined with the boundary condition (Q) , the boundary 
effects can be regarded as being formulated completely. The investigation 
below, however, tells a different story. First of all, note that a realistic 
boundary cannot be considered as having a geometrically and physically 
uniform surface. This, according to the spirit of calculus, leaves us no choice 
but to divide the boundary into many area elements, infinitesimally small in 
the theoretical sense. We are then supposed to examine how each of them 
receive and emit incident molecules. If we assume that an area element 
dS in Fig. 3 uniformly reemit the incident molecules and the distribution 
function immediately above it is known as /(r Q , v), it follows from (|6|) that, 
assuming the gas to be relatively stationary, 

/(r ,v) = /(r Pl ,v) = /(r P2 ,v). (14) 

An interesting question must arise. What is the difference between the ve- 
locity distributions at p± and at P2? Fig. 3b offers manifestation that the 
velocity distribution at p\ takes the value of f(r ,v) within the solid-angle 
range df^i while the velocity distribution at p2 takes the same value within 
dQ.2- An essential fact is that the infinitesimally small dfli is significantly 
larger than the infinitesimally small df^- It is then rather obvious that equa- 
tion (|l4|) alone cannot be deemed as a satisfactory description of the process. 
(The next section shows that introducing ^-functions into the scheme does 
not provide much help either.) 

We have seen that the path-invariance lemma, though intended to de- 
scribe the evolution of continuous distribution function, discloses many fea- 
tures that cannot be well treated in terms of continuous distribution func- 
tions. 

Before finishing this section and turning our attention to the next sub- 
ject, we wish to comment on the differences between the Boltzmann-equation 
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approach and the approach represented by formulas (H)-®, though it seems 
that formula (||) or (|6|) is "derived from" and "completely equivalent to" the 
collisionless Boltzmann equation ([!]). 

According to most textbooks, the derivation of the Boltzmann equation 
is closely related to the conservation law in the phase space //-space [|l[|. Ref. 
4 and Ref. 5, however, argue that such a conservation law is not truly 
sound. On one hand, defining fluxes through five-dimensional hypersurfaces 
(six of them can enclose a six-dimensional volume element) is an absolute 
must for the conservation law; on the other hand, such fluxes cannot be 
well defined in either physical or mathematical sense. For instance, in defin- 
ing a flux through the ordinary surface dydz we invoke that the velocity 
x = v x is perpendicular to the surface; whereas in defining a flux through 
the hypersurface dv x dv y dv z dydz we cannot invoke that the velocity x = v x 
is perpendicular to the hypersurface, because v x itself is among the five 
dimensions of the hypersurface. 

Formula (|j) can be derived independently through a Jacobian approach 
under two assumptions Q One is that the distribution function is per- 
fectly continuous at an initial time. The other is that forces acting on all 
molecules are free from dissipation (independent of speed) and free from 
fluctuation (smooth everywhere). After obtaining (||), equation (|l|) can be 
obtained by abandoning the path-information expressed by equation (||). 
Noting that the abandonment involves an information loss, we wish to say 
that the approach represented by instead of the collisionless Boltz- 

mann equation (Q), should be regarded as a more basic and more complete 
formalism. The discussion of this section, in which expressions in- 
stead of (H|), reveal the essential features of continuous distribution function, 
also substantiates the statement presented above. 

The standard kinetic theory gets heritages from the ordinary fluid me- 
chanics. It assumes that a local value of distribution function is largely influ- 
enced by its immediate neighborhood and the influence is exerted through 
two-dimensional surfaces or five-dimensional hypersurfaces. By adopting 
this assumption, the theory unanimously employs a conservation-type par- 
tial differential equation to set up its framework. Formulas (U)-(H), however, 
suggest that events of gas dynamics develop along molecular paths, which 
are one-dimensional lines in /i-space. The new picture, quite different from 
the standard one in the physical and mathematical senses, implies that even 
a gas that initially possesses a perfectly continuous distribution function will 
not behave itself like a continuous medium. We have manifestly seen this 
point in the discussion related to the quasi-discontinuity. 
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Finally, formulas offer, in a primitive manner though, different 

concepts concerning the way we approach to gas dynamics. Instead of rely- 
ing on comprehensive, but difficult and delicate, differential equations, one 
is virtually led to accepting a procedure-type approach. Such approach sup- 
posedly includes the following steps: (i) to determine paths of individual 
molecules; (ii) to formulate how a local distribution function gives a contri- 
bution to another local distribution function through a path; (iii) to integrate 
all contributions associated with all possible molecular paths. The last step 
is necessary since, according to the basic principles of statistical mechanics, 
molecular paths must involve a probabilistic nature due to the molecule- 
boundary and molecule-molecule interaction. (All points in this paper are 
given in terms of classical mechanics.) We will, by following Feynmanj^], 
call an approach of this type the path-integral approach. 

3 Evolution of discontinuous distribution function 

The last section has shown that the Boltzmann gas should be regarded 
less as a perfectly continuous medium and more as a special collection of 
individual molecules. In this section we will formulate how discontinuous 
distribution functions, each of which describes a set of molecules, develop 
along molecular paths. 

For simplicity only discontinuous distribution functions that are pro- 
duced by point-like sources and surface- like sources are of our interest. It 
is also assumed that molecules of the interested gases are free from exter- 
nal forces. If needed, the approach here can, in a straightforward way, be 
adapted to more general situations. 

In Fig. 4, a point ro is plotted at which molecules are generated (for 
whatever reasons). To make the discussion applicable later, the molecular 
emission rate p is allowed to depend on the time, velocity magnitude and 
velocity direction. The number of molecular emission can thus be expressed 
by 

p(to,r ,v,Q )dt Q dvdn , (15) 

where £lo is defined in the frame whose origin is at ro- Moving along a 
molecule's path, it is found that at a later time t the emitted molecules 
spread over the volume 

[r — ro\ 2 vdtdQo (16) 
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where r is the point in the spatial space such that 



r-r = v(t-t„). (17) 
It follows from (|l5|) and ([R]) that the density of molecules at r is 

pdvdt dn pdv 

n(r) = | 1 2 ,, , n = I jo-, (18) 

|r — r \ z v dtdil \r — r \ z v 

where dt = dt^ has been understood. Since the distribution function / at 
the point r satisfies 

/ fv 2 dvdfl = n(r) = , pdV l9 , (19) 
Jn |r-r pu 

the distribution function produced by the point-like source is 

f(t,r,v,u) = — t^t- 5{n - n ). (20) 

r — ropy 



Note that fio in ( |2C)|) takes the same direction as that of (r — ro). The form of 
( |20| ) manifests that molecules produced by a point-like source are associated 
with a discontinuous distribution function. If one tries to apply the Boltz- 
mann equation (|l[) to it, difficulties arise sharply with the differentiation 
operations. 

Now, we consider a surface-like source of molecules. It seems that after 
the surface is divided into many small area elements and the formula Eq. 
( p0|) is applied to each of them, our task is virtually accomplished. However, 
the following discussion shows that if we stick to what the standard kinetic 
theory implies, more troublesome things will emerge and no progress can be 
made. 

Let ASoi denote one of area elements on the surface and 77^0, ro, v, f2o) 
denote the molecular emission rate per unit area on ASoi- By identifying 

T](t ,r ,v,tt )AS 0i (21) 



with the molecular emission rate p in (|20D, we find the entire distribution 
function produced by the surface-like source to be 

fit, r, t,, n) = Y, 1 3 *(" - (22) 
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where i runs over all the elements on the surface. If one rewrites (22) in the 
integral form 

/(t,r,«,o) = / vd % 3 5(o-o ), (23) 

J |r — ro| v° 

the following question will be of the immediate concern. Does equation ( |23| ) 
define a normally behaved distribution function? Fig. 5 demonstrates that 
each of S7o points to a different direction while the distribution function on 
the left side of (p3| ) involves only one direction defined by O. That is to say, 
expression (p3|) cannot be integrated in the usual sense. 

It is quite interesting to look at the peculiarities that we have just en- 
countered. In dealing with our discontinuous distribution functions neither 
the usual differentiation nor the usual integration works smoothly. 

In order to find a way out of the difficult situation, we propose to use 
the following solid- angle- average distribution function such that if the exact 
distribution function f(t, r, v) is known, which exists only in a pure academic 
sense, then the average distribution function is 

f(t, r, v, AO) = -L / f(t, r, v)dfi, (24) 
AS l J ad. 

where Ail is a solid-angle range in the velocity space set by the investigator. 
In practical calculations, it is convenient to employ the spherical coordi- 
nate system of velocity, in which O is defined by the polar angle 9 and the 
azimuthal angle tfi, and the solid angle range can be expressed by 

AO ns sin0A0A<£. (25) 

On this understanding, the entire solid-angle range of velocity associated 
with a spatial point r is divided into a large number of small, but finite, 
ranges, which may or may not be equal to each other. Of course, if the 
distribution function is smooth enough and AO is sufficiently small, we may 
simply assume 

f(t,r,v,An) = f(t,r,v,n). (26) 

In the rest of this paper, we will always omit the bar notation when referring 
to such distribution function. 

In terms of practical calculations, the size of AO should be chosen prop- 
erly so that the computational work can be done efficiently and at the same 
time no significant macroscopic phenomena will be overlooked. 

As we may notice, an investigator of numerical work can take a sim- 
ilar strategy to treat distribution functions in terms of Ar and Av. The 
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difference here lies in that our solid-angle-average distribution function is 
introduced in a theoretical (analytical) consideration: the discontinuity rep- 
resented by (^) and the quasi-discontinuity represented by (|l3|) should, and 
have to, be handled under the new definition. 

It is rather important to emphasize that the distribution function intro- 
duced above seems to be an "approximate" one, but it actually represents an 
"accurate" approach in the following two senses. One is that by letting AO 
be sufficiently small, we can describe a statistical process with any desirable 
accuracy. The other is that errors related to giving up the exact distribution 
function are largely inherent to nature not to the way we approach to it. 

With help of the new average distribution function, the deterministic 
nature and the probabilistic nature of a statistical process can be kept in 
a balanced way. By using all kinds of mathematical operations, such as 
differentiation and integration, to formulate the distribution function, we 
preserve the macroscopic causality. By taking the solid-angle average, some 
of the microscopic information, in particular those related to discontinuity 
and quasi-discontinuity, are forsaken forever. 

We are now equipped to formulate the discontinuous distribution func- 
tion due to the existence of a surface-like source, as shown in Fig. 6. For 
a chosen solid-angle range AO in the velocity space there is an effective 
cone —Ail in the spatial space. A surface-like molecular source enclosed 
by the effective cone, the shaded area AS in Fig. 6, gives contributions to 

f(t,T,V,Att). 

Assuming the emission rate on the surface to be known and allowing 
expression (|23|), not normally behaved though, to represent the emitted 



molecules, we obtain from (24) 



f = -Ko / / | A n ~ (27) 

Ail J J as |r — rop v a 

By exchanging the order of the integration, we finally arrive at 



f(t,r,v, An) = — / — -— . (28) 



1 f ri(to,r Q ,v,n ) dS 
AO J as |r — ro| 2 v 3 

In the integrand of ( [28] ) ro is the position of dSo, to is equal to t — |r — tq\/v 
and Qo points to the direction of (r — ro). 

It is easy to see that the distribution function expressed by ( p8[ ) is finite 
and well behaved. 
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4 The complete path- integral formulation and its 
application 

The discussion in the last two sections has shown that the distribution func- 
tion at any specific point consists of two parts. The first part is a continuous 
one produced by the continuous distribution function existing previously; 
the second part is a discontinuous one produced by surface-like molecular 
sources. Here, we summarize the last two sections and give a comprehensive 
formulation for the dynamics of collisionless Boltzmann gas. 

In Fig. 7 we assume that the complete solid- angle- average distribution 
function f(t, r, v, ACl) is affected by a continuous distribution function and 
a surface-like source. Both the "sources" exist within the effective cone 
and constantly "emit" molecules into the region around the point r. The 
complete distribution function can formally be expressed as 

/(t,r,«,Afi) = / (0 +/ (ii)) (29) 

where fu) and Jua are produced by the aforementioned two sources respec- 
tively. 

First, we wish to determine in terms of certain specifications of ini- 
tial state. Providing the initial state of the continuous distribution function, 
denoted by f ct , is given only at one specific moment to, we have no other 
choices but to take the following approach. All relevant source points, ac- 
cording to Sec. 2, can be determined by the mapping T _1 r and by the 
effective cone — Af2. Namely, they distribute on an area element of the 
spherical surface |ro — r| = v(t — to), denoted by ASi in Fig. 7. We will refer 
to such area element as a virtual effective surface in view of that similar 
surface-like sources (boundaries) are physical ones. From (^) and (|24|), the 
continuous part of the distribution function can be written as an integral 



f (i) (t, r, v, AO) = -!- [ f ct (t, r, v, tydtt 

AO./ A< , (3Q) 

f ct (to,r ,v,Q, )U ror dS 



1 

|r - r | 2 AO J ASi 

where ro is the position of dSo, fto is in the direction (r — ro), and U ror is a 
specially defined path-clearness function such that 

_ | 1 no block along the path from ro to r . . 

r ° r ] otherwise. 
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Note that the difficulty of quasi-discontinuity discussed in Sec. 2 is no longer 
an issue: fu\ is uniquely defined and makes appropriate physical sense no 
matter how dynamically the source function f ct (to) varies on the virtual 
effective surface AS\. 

If the distribution function is known within a period of time (before t), 
we have freedom to choose the virtual effective surface. For instance, by 
choosing the surface AS[ in Fig. 7, we have, as an alternative of (| 



1 f f ct (t ,r ,v,n )U ri)r ^ 
hi) = "ao / i 12 cosadSo, 32 

where a is the angle between the direction (r— ro) and the normal of the area 
element dSo, and the time to now depends on the position of dSo, namely 
t = t - |r - r \/v. 

Then, we wish to formulate the discontinuous part of the distribution 
function. As suggested in the last two sections, boundary surfaces in a gas 
can be treated as surfacedike molecular sources, since they constantly reemit 
incident molecules. In Fig. 7, the boundary surface AS2 within the effective 
cone, called the physical effective surface, is singled out as a surfacedike 
source that can affect f(t, r, v, Afl). 

The focus is naturally on the local emission rate tj(to, tq, v, CIq) defined by 
(f2~l|). To determine it, we first consider the falling rate of incident molecules 
expressed by 

£(to,T ,Vi,£li) = f(to,r ,Vi,£li) cosOiU(cQ&Qi), (33) 

where 9i is the angle between the inward normal of the local surface and the 
incident direction of the molecules and U is a step function whose value is 
equal to unity if cos#j > and equal to zero otherwise. 

As well known, the functional relation between 77 and £ has to be ulti- 
mately measured in experiments ||, though some kinds of theoretical models 
may be of use at some stages. To see how this relation can be formulated 
empirically, consider a molecule that moves with the velocity (vi, Oj), strikes 
an area element and leaves the area element with a velocity within the range 
dvdQo, with respect to the surface element, in a certain probability. If the 
probability is denoted by 

P(vi, fli, v, Qo)dvd£lo (34) 
and its normalization takes the form 

J P(vi, Qi, v, £lo)dvd£lo = 1, (35) 
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then the emission rate of the area element, for rarefied gases, will satisfy 

i](t ,r ,v,n ) = / P(v i ,n i ,v,n )^vfdvidri i . (36) 



To verify the correctness of pq), one may integrate ( |3q ) and obtain, with 
the help of ©, 

dt dSo J r](t Q ,r , v,£l )dvdSl = dt dS J ((t ,r ,Vi,Ui)vfdvi(Kli, (37) 

which is nothing but the conservation law of the molecular number on the 
surface element. It is now clear that the purpose of such experiment should 
be to determine the functional form of P(vi, Qi,v, fo- 
under the assumption that the molecular emission rate r\ of the involved 
boundary surface has been determined (by whatever means), the interested 
discontinuous part can, by virtue of Eq. fl28|), be written as 



1 f r](t ,ro,v,n )dSo 
K ' AU Jas 2 l r - ro| « 



where the path-clearness function U TQT has been defined by (|3T 
The final result is then 

ru ao\ 1 f / Ci (*o,r ,^,^o)|cosa| 

f{t,r,v, An) = — / : U ror dS 

AttJAS[ |r — r r 

1 f r)(t , r ,v,n ) 

+ AO / H 12 3 ^r rQJ0) 

AUJas 2 \ r — ro| u 



(39) 



where the first term of the right side can be replaced by (p0[). If macroscopic 
quantities are of interest, the distribution function expressed by ( p9| ) can be 
used as a conventional one. 

At this point, it is interesting to comment on the time-irreversibility 
of the formalism given above. It is kind of well-known that the time- 
reversibility dilemma related to a dynamical process can be eliminated if 
a proper statistical average is taken. With the statistical average, two is- 
sues get involved, (i) Some pieces of microscopic information are forsaken, 
(ii) Conservative basic forces are converted into fluctuating and dissipative 
(speed-dependent) forces. Notably, the second issue mentioned above is di- 
rectly responsible for the time-irreversibility of the interested process. In this 
approach, it appears that the time-irreversibility has nicely been embedded. 
Though the pure evolution of continuous distribution function expressed by 
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(Q) or (|8|) is time-reversible, other treatments, including the definition of 
the solid-angle-average distribution function and the formulation of bound- 
ary effects, are manifestly time-irreversible. In particular, by allowing the 
molecule-boundary interaction to be determined by empirical laws (or other 
adequate statistical laws), the formulation in this paper is fully in harmony 
with the Langevin theory in which fluctuating and dissipative forces nat- 
urally arise from the interaction between a moving body and its surrounding 
molecules. 

At the end of this section, we give a brief look at what will happen 
if a calculation suggested by this paper is practically performed. Fig. 8 
offers schematic of a gas leaking out of a large container through a small 
hole. (Ref. 4 shows that the standard theory encounters many difficulties in 
treating the case.) Since the container is rather large one may assume that 
the distribution function on a surface, labeled as S in the figure, constantly 
takes the value 

/ m \ 3 / 2 / mv 2 \ fA . 

/ = "»U^) «p(-5ff). <40) 

provided that the surface S is not very close to the hole. It is obvious that 
the values of distribution function at the starting points of paths 1,2 and 3 



are all equal to the value of (40). If contributions from complex paths, such 
as path 3, are neglected, we can directly use an ordinary PC-type computer 
to calculate (|39|) on the understanding that the molecular emission rate r] 



in the formulation has been determined empirically. 



5 Conclusion 

We have set up a new formalism that exhibits many features strongly dif- 
ferent from those related to the standard theory. They can briefly be listed 
as the following. 

1. Instead of pursuing the exact distribution function, a special type of 
distribution function, called the solid-angle-average distribution func- 
tion, is introduced. The new-type distribution function can describe 
statistical phenomena with any desirable accuracy. 

2. Instead of relying on a partial differential equation, which is difficult 
and delicate, the new approach gives a set of tamable integral formulas. 
As well known, for a law of nature the integral formulation often enjoys 
advantages. 



15 



3. Instead of disregarding the path-information, the new approach bases 
its formulation on the path- information of molecules. With the path- 
information included, many sophisticated and important features of 
kinetic systems become explicable and treatable. 

4. Instead of treating the collisionless Boltzmann g completely 
continuous medium, the new approach treats the discontinuity and 
quasi-discontinuity at the very beginning. No singularity of any type 
exists with the final result of the formalism. 

5. Instead of falling into the time-reversibility paradox, the new approach 
admits its inability to formulate all microscopic details. In treating 
boundary effects, the time irreversibility has explicitly been embedded 
in the formalism. 



The complete dynamics of the Boltzmann gas, including collisional ef- 
fects, will be formulated in different works Q [11 1. Helpful discussion with 
Prof. Keying Guan is greatly appreciated. This work is partly supported 
by the fund provided by Education Ministry, PRC. 
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Figure captions 

1. Local distribution as an image of the entire space. 

2. Source points affecting the local distribution within v 2 AQAv . 

3. Local distributions near an area element. 

4. A point-like molecular source. 

5. Schematic of discontinuous distribution produced by a surface-like 
source. 

6. A surface-like molecular source within an effective cone. 

7. One physical surface and two virtual surfaces within an effective cone. 

8. Molecules leaking out of a large container. 
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Figure 1 




Figure 2 




Figure 3 




19 



Figure 4 
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Figure 5 
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Figure 7 



Figure 8 
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